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1. Introduction 

Rayleigh-Benard convection is a paradigm of a pattern forming system far from 
equilibrium. Convective fluid motion in a vessel is induced by a vertical temperature 
gradient between the bottom and top boundaries due to buoyancy forces. In dependence 
on this temperature gradient, the geometry of the experiment and the fluid properties, 
a whole zoo of instabilities has been observed ranging from laminar, spatially coherent 
convective motion over spatially ordered but temporally chaotic up to highly turbulent 
fluid motion. We refer the reader to reviews available on the topic [H [2]. 

Recently, much efforts have been devoted to the analysis of turbulent Rayleigh- 
Benard (RB) convection both by experimental as well as theoretical means [Hill]- Direct 
numerical simulations allow one to consider the dynamical and statistical properties of 
RB turbulence and the transitions between different types of flows in fine detail. 

It is obvious that the analysis of turbulent convective fluid motion has to be 
based on a combination of tools from dynamical systems theory, statistical physics, 
and the theory of stochastic processes. A necessary step is the statistical formulation 
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of the underlying basic fluid dynamic equations, which for the most simple case are 
the Oberbeck-Boussinesq equations for the velocity field u{r,t), the temperature field 
T{r,t), and the pressure field p{r,t): 

^ -^T(r,t)=Ar(r,t) 



+ U-V 



dt 

d_ 
dt 

V ■ u{r,t) 



u{r,t) 



Vp{r, t) + Ft [Au{r, t) + RaT{r, t)ej (1.1) 







The equations have been nondimensionalized using the Rayleigh number Ra = "^'^J^ , 
which is a dimensionless measure of the temperature gradient across the fluid layer 
(with thermal expansion coefficient a, gravitational acceleration g, outer temperature 
difference AT, and distance of top and bottom plate h), as well as the Prandtl number 
Pr = ^ as the ratio of kinematic viscosity u to heat conductivity n of the fluid. Thus, 
the vertical spatial coordinate obeys z G [0,1], and the boundary conditions of the 
temperature at the bottom and top plate are T{z = 0) = ^ and T{z = 1) = — |. 
For the velocity, no-slip boundary conditions u{z = 0) = u{z = 1) = are assumed. 
These equations are solved numerically by a suitably designed penalization approach, 
described in section [5j A snapshot of the temperature field is exhibited in figure [Tj 




Figure 1. Snapshot of the temperature field T{r). The green, blue and red arrows 
point in x-, y- and z-direction, respectively. The volume rendering has been done with 
the open source renderer Voreen [S]. 



The statistical analysis is based on joint probability density functions (PDFs) for 
the temperature and the velocity at a single point in space and time. The basic fluid 
dynamic equations require the validity of certain relations among these PDFs. From 
these relations, corresponding expressions relating the various moments of the fields 
can be derived. For the case of incompressible turbulence, these relations have been 
formulated by Lundgren, Monin and Novikov [HI [71 |8j, and Ulinich and Lyubimov |9] 
and are sometimes known as the Lundgren- Monin- Novikov (LMN) hierarchy. They 
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are directly related to Hopf's functional equation, which can be viewed as the basic 
statistical formulation of the Navier-Stokes equation in the Eulerian framework [10] . 
Similar relations can be derived for the corresponding Lagrangian quantities [H]. It is 
evident that an analogous treatment is feasible for the Oberbeck-Boussinesq equations. 

In the present article we will use this approach in order to analyze the single-point 
temperature probability density function for stationary turbulent RB convection. Our 
analysis combines direct numerical simulations with the relation of the LMN hierarchy 
for the single-point PDF. The result is a partial differential equation for the temperature 
PDF. The derivation of this relation is outlined in section section |2l 

The aim is to formulate an equation which characterizes this PDF. Starting from 
the Oberbeck-Boussinesq equations, we derive an evolution equation for the single- 
point joint probability density function of velocity and temperature along the lines of 
Lundgren, Monin and Novikov P, d E] for the case of incompressible turbulence. This 
equation is unclosed due to the fact that it contains unclosed expressions which can 
be related to fluid pressure, viscous dissipation and heat diffusion. However, these 
expressions can be treated by introducing conditional averages, which can be extracted 
from direct numerical simulations. This leads to a partial differential equation for the 
joint temperature-velocity PDF. A similar approach has been performed by Novikov 
[T^ [T3] . and more recently by Wilczek et al. for the PDFs of vorticity [13] and velocity 
[15] for stationary, isotropic turbulence. On the other hand, modeling of unclosed terms 
is also a possible method, as performed by e.g. Pope [IGl IT7] . As we shall indicate, 
the analysis of the evolution equation for the temperature PDF yields a comprehensive 
description of the dynamical processes in RB convection. 

The article is structured as follows: In section |2] we will derive the evolution 
equation for the temperature-velocity joint PDF. In section [3} we reduce the joint 
PDF to the temperature PDF and make use of statistical symmetries to cut down 
the complexity of the evolution equation. Then we present a descriptive way to deal 
with this equation involving the method of characteristics. Section |4] briefly discusses 
connections to the Nusselt number, relevant for the heat transport. These theoretical 
results are complemented by results from direct numerical simulations, which will be 
discussed in section |5| followed by a summary in section [6j 

2. Single- Point PDF 

We are interested in the joint temperature-velocity probability distribution f{T,v;r,t) 
and want to derive the corresponding evolution equation. Formally, the probability 
density function is obtained as a suitable average over the so-called fine-grained 
probability distribution 

/(r, v; r, t) = 6{r - T{r, t)) 6{v - u{r, t)). (2.1) 

It is important to distinguish between sample space variables r, v and the corresponding 
realizations of the temperature and velocity fields T{r,t), u{r,t). Therefore, one could 
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think of / as the PDF of one particular reahzation of the fields. Also, the notation of the 
arguments in f{T,v;r,t) emphasizes the difference between the sample space variables 
r, V and the coordinates r, t - the PDF is normalized with respect to the sample space 
variables, the coordinates are just parameters. 

The full probability density function is now obtained as an ensemble average over 
all possible realizations of the temperature and velocity fields: 

/(r, v; r, t) = {/{r, v; r, t)) = {6{r - T{r, t)) 6{v - u{r, t))) (2.2) 

The brackets (■) denote the ensemble average, in contrast to spatial averages (•)y and 
{■)a over the whole fluid volume, or a horizontal plane at height z, respectively. 

The definition of the fine-grained PDF (2.1) can be differentiated with respect to 
the space and time variables, giving 



_d_ 



d dT{r,t) 

Or dxi 



/ + V„ 



du{r, t) 

dxi 



f 



(2.3) 



Note 



as the spatial gradient, and an analogous equation for the te mpor al derivative 
that the operators and V,; act on /. Now, multiplying (2.3) by Ui and adding the 
temporal derivative allows us to make use of the basic Oberbeck-Boussinesq equations 
(1.1 ). This results in the desired evolution equation for the fine-grained PDF: 

/ du 



d_ 
dt 



f + u-Vf=- 



d_ 
'd_ 



dT 



\dt 

Pr (Aw + RaTe^)]/ 



+ u ■ VTj f + V^-{— + u-Vu]f 



(2.4) 



dt 

-ATf- V„ 

Performing the ensemble average of this equation in order to arrive at an equation 
for the full PDF (2.2), one encounters the closure problem of turbulence, since the 
unclosed averages (AT/), (— Vp/) and (Auf) show up. The LMN hierarchy ansatz 
would mean to treat these terms via a coupling to the two-point PDF; the evolution 
equation of this two-point PDF would in turn introduce a coupling to the three-point 
PDF, and so on 

Instead of introducing this hierarchy of coupled evolution equations for the multi- 
point PDFs, our strategy [HI [T5 ] is to express the unclosed terms as conditional 
averages, since these are accessible to direct numerical simulations. The result is a 
partial differential equation governing the joint temperature-velocity PDF and relating 
its shape as a function of space point r to the conditional averages. The functional form 
of these conditional averages is a signature of the underlying dynamical processes of RB 
convection. 

Introducing the conditional averages 

(AT/) ={AT\r,v,r,t)f 

{-Vpf) = {-Vp\T,v,r,t)f (2.5) 
(Auf) = {Au\T,v,r,t)f 

we arrive at the following relation for the single-point probability distribution 
f{r,v;r,t): 
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-f + vVf= ~—{AT\T,v,r,t)f 

- V„ ■ [(- Vp|r, V, r, t) + Ft ((Aw|r, v, r, t) + Ra re,)] / (2.6) 

Since at the boundaries of the RB cell the velocity and temperature fields are 
statistically sharp quantities, the probability distribution has to obey the conditions 
/(r, V] z = = 5 (t — \) 6{v) and /(r, v; z = 1) = 5 + ^) 5{v) for arbitrary x,y. 

We note in passing that a different version of the evolution equation can be obtained 
by introducing the Laplacian of the PDF. Specializing in the case Pr = 1, it is possible 
to re-express the conditional averages {AT\T,v,r,t) and {Au\T,v,r,t) via the relation 

A/= -^{ATf)-V.-{Auf) 

or 

+ |^((VT)V1 + 2^(Vr . VuJ) + ^(V«. . VuJ). (2.7) 

Here and in the following, Einstein's summation convention over repeated indices is 
used. Again, we can introduce conditional expectations, where the arguments of the 
conditional averages have been abbreviated as ic = T,v,r,t. The resulting evolution 
equation reads 

^/ + „ ■ V/ = A/ - V„ ■ |(-Vp|*) + Ra.TeJ / - g^{(Vr)»/ 

A somewhat more complicated equation holds for Pr 7^ 1. This relationship shows that 
the single-point joint temperature-velocity PDF f{T,v;r,t) is essentially determined 
by the conditionally averaged dissipation-like terms ((VT)^|7^r), (VT-Vmj|*) and 
{Vui ■ Vujl-k) as well as the conditional pressure gradient (— Vp|^). 

3. Single-Point Temperature PDF and Implications of Statistical 
Symmetries 

In the following we shall restrict our attention to the reduced temperature probability 
distribution and its evolution equation. As it will turn out, already this equation gives 
insight into the connection between RB dynamics and temperature statistics, besides 
obviously describing the temperature statistics itself. Also, because the final evolution 
equation of temperature PDF involves scalar valued functions only, a numerical approach 
is easily feasible. 

The reduced temperature PDF is obtained by integrating out the velocity part: 

h{r;r,t) = J d\ f{T,v;r,t) (3.1) 



Starting from equation (2.6) we obtain the simple equation 



Q-h + V- {u\t, r,t)h = - — {AT\t, r, t)h, (3.2) 
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where we have performed the integration with respect to the velocity v. Thereby, we 
had to introduce the conditional averages {u\T,r,t) and (AT|r, r,t). Alternatively, one 
could re-enact the derivation performed in section [2] for the fine-grained PDF of the 
temperature /i(r; r, t) = 6{t — T{r, t)). 



Analogous to (2.8), we can derive a further relation using the identity 

d - 92 



Ah = {Ah) = - — {ATir,t)h) + —{{VT{r,t)fh). (3.3) 



This relation also follows from (2.7). With this equation, we can now summarize the 
equation for the temperature PDF /i(r; r,t) in the form 



a, _ 92 



^h + W ■ {u\T,r,t)h = Ah - ^{iyTf\T,r,t)h. (3.4) 

Here, the conditional average of the term (VT)^ comes up, which is related to the 
Nusselt number. Details will be discussed in section HI 

The evaluation of the conditional averages appearing in (3.2) is greatly simplified by 
considering convection that is statistically stationary in time and has periodic horizontal 
boundaries (i.e. is homogeneous with respect to the horizontal coordinates). Under the 
assumption of a statistically stationary flow, the PDFs and therefore the conditional 
averages cannot depend on the time variable; also the dependency on the horizontal 
coordinates drops out. So instead of dealing with statistical quantities that depend on 
r and t, we simply have to retain the z-dependence. 

Let us first consider the determining equation for the temperature PDF h[T\ z) in 
the form 

^^{u,\T,z)h{T-z) + ^{AT\T,z)h{T-z) = Q, (3.5) 



obtained from (3.2). This equation in principle has to be solved together with the 
boundary conditions h{T] 0) = ^{j ~ \^ /i(r; 1) = (5^r + and with appropriately 
modeled expressions for the conditional averages {uz\t,z) and (AT|r, z). This direct 
approach will not be conducted in the present paper though, but can be taken as a 
starting point for future modeling. 

Instead, this first order partial differential equation (PDE) can be analyzed with 
the help of the method of characteristics [18]. Applying this method, one can find 
curves {t{s),z{s)^ in the r-z-phase space parameterized by s along which the PDE 
(3.5) transforms into an ordinary differential equation which can be integrated. This 
approach will be sketched in the following. 

Writing h{s) = /i(r(s); z(s) ) and calculating the derivative of it gives 



d / , , , ^\ dhdr dhdz , , 

-/,(rW,.w) = + (3.6) 

The PDE ( 3.5[ ) is re-expressed in the form 

{AT\t,z)^ + {u.\t,z)^^ = - (^{AT\t,z) + §-^{u.\t,z)\ h. (3.7) 
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Comparing these two equations identifies the characteristic curves as solutions of 
d 

d^ 
d 

d^' 



-t={AT\t,z) 
-z = {uz\t,z) 



(3.8) 



Along these curves, the PDE (3.5) becomes 

d , f d , , ^, , d , , , , , 



which can be integrated to 

h{s) = h^so) exp 



ds{ ^{AT\T,z) + -^{u,\r,z) 



so 



(3.9) 



(3.10) 



This equation describes the evolution of the PDF along a trajectory ( r(s), z{s) ) starting 



at point yr^so), z{so) j in phase space. A particularly appealing property of this 
formalism is that it allows to interpret the statistical results in an illustrative manner, 
because the characteristics, i.e. trajectories in r-z-phase space, show the evolution of 
the "averaged" physical process. 

It is tempting to interpret the characteristics as a kind of Lagrangian dynamics 
of a tracer particle inside the RB cell. However, the dynamic of a tracer particle is 



stochastic, whereas the characteristics defined by (3.8) describe purely deterministic 



trajectories and, thus, take the stochastic properties into account only in an averaged 
way. In a sense, the characteristics describe the averaged evolution of an ensemble of 
fluid particles that are defined by their initial condition in the r-2;-plane. 

Thinking of turbulent RB convection with some physical intuition, one can 
expect certain features from the statistical quantities introduced in this section. The 
conditionally averaged vertical velocity {uz\t,z) should show positive correlation with 
the temperature, i.e. it should mirror the well-known fact that hot fluid rises up and 
cold fluid sinks down. Also the no-slip boundaries should be recognizable for z ^ and 
z ^ 1, respectively. The absolute value of the heat diffusive term {AT\t,z) should be 
highest near the boundaries because of the sharp change of the temperature profile. As 
the characteristics in a way describe the average path a fluid particle takes through t-z- 
phase space, the typical Rayleigh-Benard cycle of fluid heating up at the bottom, rising 
up, cooling down at the top and sinking down again should find its correspondence in 
the statistical quantities describing the evolution of the PDF. Actual numerical studies 
of these quantities will be discussed in detail in section [5j 

In the 1990s, an approach similar to the presented one was undertaken by Yakhot 
and followers, where PDFs of various quantities in a stationary flow were expressed as 
integrals over conditionally averaged variables. Considered quantities included passive 
scalars by Sinai and Yakhot [19], active scalars, such as temperature fluctuations by 
Yakhot [20] and temperature increments by Ching [21], and even general functions of 
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an arbitrary quantity measured in the flow by Pope and Ching [221123]. The conditional 
averages included spatial and temporal derivatives of these quantities. 

In contrast to these works we do not assume homogeneity in the spatial coordinates, 
i.e. we still have the z-dependency present in our PDF equations. This allows us to 
discuss the PDFs with respect to the ^-coordinate and observe qualitatively different 
statistics in different regions of the flow. As a result, we are able to see the differing 
behaviour of bulk and boundary parts of the convection cell and how these are connected 
to the conditional averages and their dependency on the vertical position. This has not 
been addressed in the literature. 

4. Connection to the Heat Transport 

The Nusselt number Nu as the ratio of convective to conductive heat transfer plays a 
key role in the analysis of RB convection. It serves as a measure of how efficient heat 
can be transported through the convection cell. 

Though we did not consider the Nusselt number so far, we note that there is 
an interesting connection between Nu and the conditional averages that appear in 
our calculations. Because the Nusselt number may be defined as the volume average 
Nu = ((VT)^)y in non-dimensional units, the following expression comes up: 

Nu = ^ j d^r dr {{VTf\T,r)h{r;r) (4.1) 

Therefore, ((VT) |r,r) can be viewed as a conditional Nusselt number density. We 
point out that this quantity should be of considerable interest for the evaluation of 
theories concerning the Rayleigh number dependency of the Nusselt number based on 
a decomposition of the heat transport into bulk and boundary contributions, which 
underlies the Grossmann-Lohse theory [21] outlined in the review of Ahlers et al. [1]. 
Also, this term is linked with the temperature dissipation rate, which is discussed in 
[25] with respect to temperature PDFs. 

In a similar manner, we can employ the temperature-velocity joint PDF /(r, v; r) 
to derive an equation relating Rayleigh and Nusselt number. From the relation 
Ra(Nu— 1) = {(Vu)'^)y it is straightforward to see that 

jd^rdTd\f{T,v;r)[Ra{{{VTf\T,v,r)-l)-{{Vuf\T,v,r)] =0. (4.2) 

These two exact relations underline the importance of the conditional averages of (VT)^ 
and (Vu)^ that naturally come up in our derivations. 

5. Numerical Results 



The benefit of our theoretical approach is that we can easily provide it with 
measurements and data in the form of numerical results. To this end, we solve the basic 
Oberbeck-Boussinesq equations (1.1) with a standard dealiased pseudospectral code on 
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a three-dimensional equidistant Cartesian grid with periodic boundary conditions. For 
an introduction into this topic the reader is referred to |26| [27]. 

Periodic boundaries are required in horizontal direction, but in vertical direction 
Dirichlet conditions for velocity and temperature (i.e. no-slip boundaries of constant 
temperature) are needed. They are enforced by a volume penalization ansatz |2S t l2 ^ IBD]: 
The fluid domain VL = [0, L^] x [0, Ly] x [0, 1] C is embedded in a computational 
domain that is extended by a layer of thickness d in z-direction, = [0, L^] x [0, Ly] x 
[—d, 1+d]. Inside the fluid domain Q G Qc the unaltered Oberbeck-Boussinesq equations 
are solved, while in the appended extra regions Qc\^ a strong exponential damping 
{—-u and —-6, respectively, with ^ 1) is added to the evolution equations (1.1) of 
velocity and temperature that damps the fields to zero. By simulating the deviation from 
the linear temperature profile, 0{r,t) := T{r,t) + {z — ■^), instead of the temperature 
itself, the desired boundary conditions read u = and ^ = for z = and z = 1. This 
change of variables T ^ 9 allows us to make use of the volume penalization approach 
in a straightforward manner. 

The reason we choose this numerical scheme instead of often used Chebyshev-based 
codes, described for example in and used in e.g. pij [52] . is because it allows for 
almost arbitrary shaped boundaries and sidewalls. Although this feature is not used in 
the present paper due to the required horizontal homogeneity, it even allows to simulate 
cylindrical vessels on a (numerically cheap) Cartesian grid. A more detailed report on 
this will be published in the future. 

Our theoretical derivation relies on the concept of ensemble averages. Of course, 
through our numerics we can only access a finite subset of all possible ensemble members. 
So due to the statistical symmetries and by assuming ergodicity, the ensemble average is 
substituted for a combined volume-time average in the numerics. Likewise, the volume- 
averages and introduced earlier are actually evaluated as combined volume and 
time averages. Here, time averaging means averaging over 1250 statistically independent 
snapshots of the fields. The line plots below only show the parts of the statistical 
quantities where the statistics converged, i.e. where a significant number of events was 
obtained. 

The simulation was conducted for the parameters Ra = 4.33 x 10 and Pr = 1. The 
computational domain Oc is resolved with x Ny x Nz = 256 x 256 x 192 gridpoints 
on an equidistant Cartesian grid, where the fluid domain i7 C fic is represented by 
256 X 256 X 128 gridpoints. Thus, the aspect ratio is F = 2, with the two horizontal 
dimensions being identical. The Nusselt number is estimated to Nu = 24.2. We mention 
that we repeated the analysis presented below with a simulation of aspect ratio F = 4 
and found basically the same features, due to the fact that for sufficiently large F the 
periodic boundaries do not have a significant infiuence on the fiow. 

Figure [2|^a) shows the color-coded temperature PDF h{T] z) in a logarithmic plot, 
where the black and gray lines indicate the mean temperature profile {T{z))a and 

the square root of the centralized second moment, ( (T — {T{z))a) ) , i.e. standard 
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Figure 2. (a) The mean temperature {T{z))a and a color plot of the logarithm of the 
temperature PDF /i(r; z) as a function of z. The additional solid gray lines mark the 
contour line for {T{z))a ^ {{T ~ {T{z)) a)^)1['^ and indicate the standard deviation of 
temperature at height z. The horizontal dashed lines indicate the positions of slices in 
T-direction at fixed height and are located at z £ j, 4(5t, 2St, St}, where St — 
is the thermal boundary layer thickness, (b) The logarithm of the temperature PDF 
/i(r; z) for different values of z. The upper abscissa is scaled in units of the globally 
taken standard deviation of temperature, T^ms = a/ {T'^)v- 



deviation. Figure [2]^b) shows slices in r-direction indicated by the dashed hnes in 
figure [2]^a). In the color plot, one clearly observes the sharp change of the temperature 
PDF from a 5-function at the boundaries across the boundary layer to a shape exhibiting 
larger tails in the bulk. In addition to these tails, another feature of the PDF is the 
hump close to the r = 0-line. This hump corresponds to the most probable value 
of the temperature. One expects two different dynamical features to be responsible 
for this special shape; a tempting explanation would be to attribute the hump to the 
background temperature field of mean temperature, and account the wings for large |r| 
for plumes that carry fiuid that is much colder or hotter than the surrounding fiuid. An 
evidence for this is the shape of the PDF close to the bottom or top boundaries (but 
still outside the boundary layers): hi z = 4(5t, though the most probable temperature 
value is moved slightly towards lower temperatures, the PDF exhibits a large tail at high 
temperatures. The interpretation is that mostly cold fiuid gathers in the lower regions 
of the bulk, being almost at rest (compare the region of (u^lr, z = 45^) in figure |3](b) 
corresponding to the hump), while very hot fiuid is a more rare event, because hot fiuid 
is convected away quickly due to plume dynamics. The reason why the hot fiuid takes 
greater temperature values than the cold fiuid (in terms of absolute value) is that very 
cold fiuid detaching from the top plate already heats up on its way down. 

PDFs of the same shape for Rayleigh numbers of the same order are reported 
in j25], where also the dependence on the vertical coordinate is taken into account. 
The experimental data in [33| [2T] shows a more pronounced exponential shape of the 
temperature PDF, which can be attributed to the difference in the Rayleigh numbers 
which are several orders of magnitude above ours; the numerical data in [25] suggests 
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Figure 3. (a) Color plot of the conditional velocity field {uz\t,z), with dashed lines 
as in figure [2| Due to the non-dimensionalization, the velocity is given in units of the 
heat diffusion velocity, i.e. the velocity with which heat would be transported from 
plate to plate by pure heat diffusion, (b) The conditional velocity field {uz\t,z) as a 
function of r for various values of z. 
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Figure 4. (a) Color plot of the conditional heat diffusion field (AT|t, z), with dashed 
lines as in figure [2] (b) The conditional heat diffusion field (AT|r, z) as a function of 
T for various values of z. 



that the PDFs become more exponential with increasing Rayleigh number. 

Figure [3] and |4] exhibit the conditional averages introduced in section |3} One can 
clearly observe the features that were suggested in the aforementioned section. The 
conditional vertical velocity {uz\t,z) is high (low) for hot (cold) fluid respectively, and 
the no-slip boundary conditions manifest in the fact that {Uz\t,z) is close to zero for 
z ^ and z ^ 1. Additionally, one observes a stripe close to the r = 0-line of almost 
vanishing vertical velocity which coincides with the reddish core (the hump, i.e. the 
most probable value) of the temperature PDF in figure [2]^a). The interpretation is that 
fluid that is as hot as the mean temperature is neutrally buoyant and neither moves up 
nor down. Another striking feature is the sudden increase of the vertical velocity for 
high r near the boundary layer, i.e. for z = 6t, which we attribute to rising plumes that 
detach from the hot bottom plate. Again, it must be stressed that these interpretations 
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hold in an averaged sense. 

Figure |4] shows that the conditional heat diffusion term (AT|r, z) is (in terms of 
absolute value) highest at the boundaries, with the term being positive (negative) at the 
hot bottom (cold top) plate. On the contrary, in the bulk the absolute value is high (low) 
for very cold (hot) fluid, i.e. in the wings of the temperature PDF. Additionally, the 
r-slice near the boundary in figure |4](b) shows an under- and overshoot. The connection 
of these unique features to the RB dynamics has yet to be understood. 




-1/2 1/2 



Figure 5. Color plot of the temperature PDF /i(t; z) together with the vector field 
identifying the characteristics. 



By combining the two aforementioned conditional averages to the vector field (3.8) 
that defines the characteristics as suggested in section |3| one arrives at the vector field 
depicted in figure [5] - one of our central results. It is easy to interpret this graph 
by tracing the vector field; one can qualitatively reconstruct the typical RB cycle of 
fluid heating up at the bottom, rising up while starting to cool down, cooling down 
drastically at the top plate, falling down towards the bottom plate while warming up 
a bit and heating up again at the bottom. It is especially illustrative to see that the 
main contribution of cooling and heating (i.e. biggest movement in r-direction of phase 
space) takes place near the boundaries, highlighting the importance of the boundary 
layers, while obviously the biggest movement in 2;-direction occurs in the bulk. 

Yet one has to consider, for example, that although hot fluid rises up very quickly 
(referring to the vectors pointing upwards at the right side in figure [s]), this does not 
contribute much to heat transport because these events occur rarely, as indicated by 
the temperature PDF shown along with the vector field governing the characteristics. 

In figure pi the conditional heat dissipation rate ((VT)^|r, 2) is shown, which can 



be interpreted as a Nusselt number density according to (4.1). The noteable features 
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Figure 6. (a) Color plot of the conditional field ((Vr)^|T, z), with dashed lines as 
in figure [2] The color axis is scaled logarithmically, (b) The field ((Vr)^|T, z) as a 
function of r for various values of z. 

are a pronounced minimum near the most probable temperature, again coinciding with 
the reddish core of the PDF. Also, in the boundary layer this quantity attains huge 
values (note the logarithmic scaling in the color plot), which highlights the fact that the 
boundary layer contributes much to the heat transport. A similar shape of a related 
quantity is reported in [25j (using the deviation from the mean temperature profile 
instead of the temperature itself), though there the conditional average is taken over 
the whole fluid volume and is hence lacking the z-dependency. 



6. Summary 

In the present work, we have analyzed the single-point temperature PDF on the basis 
of the Lundgren-Monin-Novikov hierarchy by truncating the hierarchy on the first level 
via the introduction of conditional averages. 

We have first derived the evolution equation of the full joint PDF of temperature 
and velocity. Then we focused on the temperature PDF only, which is the central point 
of our paper, and obtained an evolution equation for it by reducing the joint PDF 
equation. We assumed rather weak symmetry conditions of statistical stationarity in 
time and homogeneity in lateral spatial directions; these conditions should be fulfilled 
at reasonably high aspect ratios even for closed vessels, i.e. are a good approximation 
of experimental setups in the bulk of the flow. Under these symmetry considerations, 
the evolution equation of the temperature PDF becomes fairly simple. The arising 
conditional averages of temperature diffusion (AT|r, and vertical velocity {uz\t,z) 
are estimated by direct numerical simulations using a suitably designed penalization 
approach, and features of them are discussed. It shows that expected features such as 
properties of the temperature and velocity boundary layers, correlation of temperature 
and velocity and so on are related to the form of these conditional averages that naturally 
come up in our derivations. 

The evolution equation of the temperature PDF is readily treated by the method 
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of characteristics. Due to the apphed symmetry conditions, the phase space which 
describes our system becomes two-dimensional, spanned by temperature r and vertical 
coordinate z. Because of this reduced dimensionality of the system, the method of 
characteristics yields a descriptive view of the RB dynamics, resulting in the vector 
field describing the evolution in r-2;-phase space. The characteristics, i.e. trajectories 
in r-2-phase space, are found to reproduce the typical cycle of a fluid parcel. The 
regions of the main transport in r-direction have been identified as the boundary layers, 
while the major movement in z-direction takes place in the bulk. This highlights the 
importance of the boundary layers to the heat transport. The relation of heat transport 
in terms of the Nusselt number to the conditional averages introduced in our derivation 
is briefly discussed, leading us to the definition of a Nusselt number density. It would 
be very interesting to obtain the statistical quantities describing the evolution of the 
PDF directly from experiments, e.g. from measurements of instrumented particles as 
described in [M] . 

Future efforts will be to not only use the characteristics as an illustrative way 
to describe the mean movement in phase space, but actually calculate the PDF of 
temperature from the integral representation (3.10). Also, modeling of the conditional 
averages, which are up to now estimated from direct numerical simulations, might be 
feasible; an intermediate step would be to discuss the quantities not in the turbulent case, 
but close above the bifurcation from heat conduction to convection, where analytical 
solutions of temperature and velocity fields are available. Though an easy illustration 
in the form of trajectories in two-dimensional space will not be achievable in the case 
of the joint PDF, this approach is nevertheless promising and planned for the future, 
because we hope that already the form of the conditional averages will give insight into 
the connection of the statistics to the RB dynamics. An intermediate step would be to 
concentrate on the joint PDF of temperature and vertical velocity, which should among 
others relate to the dynamics of plumes. 
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